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The strong coupling constant as is determined using recently re-analysed e + e~ annihilation 
data collected by the JADE experiment at y/s =14 to 44 GeV. The measurements are based 
on C(q|)+NLLA predictions for various event shape observables. The calculations are found 
to describe reliably data at the lowest energies of the e + e~ continuum where non-perturbative 
contributions become important. The results for as are in good agreement with the QCD 
expectation for the running of the strong coupling constant. This is the first determination of 
as at ^/s = 14 and 22 GeV based on resummed QCD predictions. 



1 Introduction 

Tests of Quantum Chromodynamics (QCD) substantially benefit from e + e~ annihilation ex- 
periments at lower centre-of-mass energies yfs since the characteristic energy evolution of the 
theory is expected to become moremanifest towards decreasing y/s. The re-analysis of data 
collected by the JADE experiment Em, as counterpart to the LEP data, has been shown to be a 
valuable effort. Recently, data at energies down to y 7 ? = 14 GeV could be employed in state- 
of-the-art QCD studies due to the successful resurrection of the original JADE simulation and 
event reconstruction software. 

In the last decade, significant progress has been made in the theoretical calculations of event 
shape observables serving as powerful tools to investigate perturbative and non-perturbative 
aspects of QCD particularly at PETRA energies. This analysis focuses on as determinations 
based on resummed calculations for event shapes which are applied here for the first time at 
y/s = 14 and 22 GeV. Event shape data have also been used to assess the performance of various 
QCD event generators tuned to LEP data. 



2 Data Samples and MC Simulation 

The JADE detector H was operated from 
1979 to 1986 at the PETRA e+e" collider. 
It was designed as a hybrid 4-7r-detector 1 
to measure charged and neutral particles. 
The studies presented here are based on 
multihadronic data samples (with number w 
of events in brackets) at = 14.0 (1734), 
22.0 (1390), 34.6 (14372, data taken 1981- 
82), 35.0 (20688, data taken 1986), 38.3 
(1587) and 43.8 GeV (3940). Simulated 1 
data were generated using the QCD event 
generators Pythia 5.7, Ariadne 4.08 
and Herwig 5.9 a combined with the io 
JADE detector simulation. We adopted 
the parameter sets used by the OPAL ex- 
periments to describe e + e~ data at y/s = 
M z o . We also considered a predecessor ver- 
sion Jetset 6.3 H used in former JADE 
studies since it was shown to describe 
e + e~ hadronic final states. Comparisons of io 
the simulated and measured distributions 
of various integral and spectral quantities 
generally gave a good description of the 
measured data. So the simulation can be 
used to correct for detector effects. 
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3 Event Shapes at PETRA Energies 
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Figure 1: Hadron level distributions for C at y/s — 14 to 
44 GeV compared with the predictions of various QCD event 
generators (left). The error bars denote the total errors. 
Also shown are the detector correction factors Ki and the 
normalised differences Si between model and data (right). 



From the data passing the multihadronic 
selection criteria!, the distributions of the 
event shape observables thrust 1 — T, heavy 
jet mass Mjj, total and wide jet broadening 
Bt and Byy, C parameter and the differ- 
ential 2- jet rate V23 based on the Durham 

scheme are calculated (cf. El). The distributions are corrected for the limited acceptance and 
resolution of the detector and for initial state photon radiation effects using Pythia for the 
standard correction. Since mass effects due to the electroweak decay of heavy b-hadrons faking 
gluon activity in the 3-jet region are crucial at \fs = 14 and 22 GeV, we take the contribution 
e + e~ — > bb as an additional background to be subtracted from the distributions. 

As an example, the resulting hadron level data distributions for C are represented by Fig. [I]. 
For comparison, the respective distributions predicted by the QCD models based on u-, d-, s- and 
c-flavoured events are shown. In case of Pythia, there is generally a good agreement between 
the data and the model over the whole kinematic range of the observables. The performance of 
Ariadne and Herwig is more moderate at 14 GeV and improves at increasing c.m.s. energies. 
Herwig significantly underestimates the peak region of the distributions. In contrast, the 
JADE-based Jetset version fits the 14 GeV data but increasingly deviates from the data at 
higher energies. The prediction of CojETsi is clearly disfavoured by the lower energy data and 
remains worse also at higher energies. Thus we do not consider this model for our as studies. 
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Figure 2: Fits of C(a|)+NLLA to Bt at y/s = 14 to 44 GeV 
(left). Hadronisation corrections and model uncertainties 
are shown as well as the dependence of the results for as and 
the x 2 /d.o.f. on the variation of the fit range (right). 
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Figure 3: Results for qs( v / s) at = 14 to 44 GeV 
derived from 0(a|)+NLLA fits to 1 - T, M H , Bt, Bw, 
C and 7/23- The inner error bars denote experimental 
uncertainties, the outer error bars are the total errors. 



4 Measurements of as 



The determination of as is based on a combination of an exact QCD matrix element calculation 
O(ag)0 intended to describe the 3-jet region of phase space and a next-to-leading- logarithmic 
approximation % (NLLA) valid in the 2-jet region where multiple radiation of soft and collinear 
gluons from a system of two hard back-to-back partons dominate. We perform x 2 -hts of the 
theoretical predictions corrected for hadronisation effects. For the main results, we use the ln(i?)- 
matchingB for the perturbative prediction with the renormalisation scale factor = \xj 1 yfs = 1 
and Pythia for the estimation of non-perturbative contributions. As an example, Fig. |2| shows 
the fitted predictions for Bt- We generally observe stable fits and good agreement with the data 
at all c.m.s. energies with x 2 /d.o.f. ranging from about 0.2 to 2.0. In case of B\y, a significant 
excess of the theory over the data in the 3-jet region of the distributions is present. 

In principle we follow the procedure inB to estimate experimental and theoretical systematic 
errors but include additional MC modelling uncertainties. Hadronisation effects and uncertain- 
ties increase significantly at 14 GeV. Experimental errors are under control for all data samples. 
On the basis of fit and experimental errors, the individual results are consistent with each other 
within 1-2 standard deviations. For each c.m.s. energy, the as-results of the six observables were 
combined using the weighted mean methodlil. The final results obtained are listed in Tab. |]. The 
total errors are dominated by higher order uncertainties. At ^/s = 14 and 22 GeV, hadronisation 
uncertainties are very large but still of the same order as the renormalisation scale uncertainties. 



5 Conclusions 

Resummed QCD theory fits event shape data 
measured with JADE well down to = 
14 GeV and allow consistent determinations of 
as • The values obtained at the lowest energies 
are affected by large hadronisation uncertain- 
ties. The LEP-tuned Pythia Monte Carlo 
used for the estimation of non-perturbative 
effects is surprisingly well capable of describ- 
ing many aspects of e + e~ hadronic final states 
at PETRA energies. The as results obtained 
here as well as in similar analyses at higher en- 
ergies based on resummed event shapes (Fig. 
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Figure 4: as-values derived from this analysis compared 
with corresponding results at higher energies. Also shown 
is the QCD expectation for the world average of ast£l. 



QCD expectation for 
Ay 2 fit of the 0(a 3 s ) 
prediction to these values taking only experi- 
mental errors into account yields as(M z o)= 0.1213 ± 0.0006 with x 2 /d.o.f.=8.3/ll. Even con- 
sidering total errors, the hypothesis of a constant value of as is disfavoured by a fit probability 
of ps 10~ 5 . The JADE data significantly improve the verification of QCD on basis on e + e data. 



Table 1: Preliminary results for as derived from the individual results using the weighted average method. 



[GeV] 




lit exp. 


hadr. 


higher ord. 


total 


14.0 
22.0 

34.6 ('82) 
35.0 ('86) 
38.3 
43.8 


0.1704 
0.1513 
0.1409 
0.1457 
0.1397 
0.1306 


±0.0051 

±0.0043 
±0.0012 ±0.0017 
±0.0011 ±0.0020 
±0.0031 ±0.0026 
±0.0019 ±0.0032 


+0.0141 
-0.0136 

±0.0101 
±0.0071 
±0.0076 
±0.0054 
±0.0056 


+0.0143 
-0.0091 
+0.0101 
-0.0065 
+0.0086 
-0.0057 
+0.0096 
-0.0064 
+0.0084 
-0.0056 
+0.0068 
-0.0044 


+0.0206 
-0.0171 
+0.0144 
-0.0121 
+0.0114 
-0.0093 
+0.0125 
-0.0101 
+0.0108 
-0.0087 
+0.0096 
-0.0080 
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